Large-scale lognormality in turbulence modeled by Ornstein-Uhlenbeck process 
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Lognormality was found experimentally for coarse-grained squared turbulence velocity and veloc- 
ity increment when the coarsening scale is comparable to the correlation scale ol the velocity (Mouri 
et al. Phys. Fluids 21, 065107, 2009). We investigate this large-scale lognormality by using a simple 
stochastic process with correlation, the Ornstein-Uhlenbeck (OU) process. It is shown that the OU 
process has a similar large-scale lognormality, which is studied numerically and analytically. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The lognormal distribution appears in a wide range of 
natural and social phenomena (see, e.g., [l|). In fluid tur- 
bulence, it is well-known that the distribution is an im- 
portant consequence of the Kolmogorov's 1962 theory for 
modeling fluctuations of the energy cascade rate across 
scales 0| • In this 1962 theory, we have a very clear picture 
why the lognormal distribution was the first candidate of 
the cascade fluctuations in his refined phenomenology. 
Namely, the energy cascade at high Reynolds number 
can be modeled as a multiplicative process consisting of a 
large number of independently and identically distributed 
random variables. This large number is important for the 
central limit theorem to be applicable to the logarithm 
of the multiplicand. 

There is a different example of lognormally distributed 
variables in turbulence. Laboratory experiments of tur- 
bulent boundary layers in 1980's suggest that spanwise 
separations between the low-speed streaks follow the log- 
normal distribution [H, H|. In this case, the underlying 
mechanism of the lognormally distributed streaks is not 
as clear as in the Kolmogorov 1962 phenomenology be- 
cause a multiplicative structure for the streaks is not 
found immediately. 

In this paper we consider yet another lognormal ex- 
ample in turbulence, which was recently found experi- 
mentally [H, Q . In Ref. Q , it is observed that the coarse- 
grained squared velocity and squared velocity increments 
over distance r 



Ur(xc) = J^J 
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are lognormally distributed when the averaging scale R 
is set to R ~ L a . Here L u is the correlation length 
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defined as the integral scale of the velocity correlation 
function (u(x + r)u(x)} and the velocity u can be either 
longitudinal or transverse velocity component. Since the 
lognormal property holds when R is comparable to the 
large scale L u , this property is called "large-scale lognor- 
mality" . For the asymptotic regime i? 3> L u , as a usual 
result of the central limit theorem, the distribution of the 
coarse-grained quantities become closer and closer to nor- 
mal (Gaussian) distribution. Notice that the large-scale 
lognormality, which we consider here, concerns the range 
of the coarse-graining scale R ~ L u , which is different 
from the final Gaussian state in R 3> L u . 

The large-scale lognormality was observed in the ex- 
periments of grid turbulence, turbulent boundary layers, 
and turbulent jets, suggesting its universality [6J. But 
again as in the case of the streaks in the boundary layer 
turbulence, why the coarse-grained turbulence data are 
lognormally distributed is not clear. Our final goal is to 
understand its mechanism why they are so. 

As a first step to this goal, we address the follow- 
ing question: Is the large-scale lognormality a general 
property of correlated random variables, not restricted 
to large-scale fluctuations of turbulence? The previous 
results Q lead us to recognize that this balance between 
R and L is most important for the lognormality. In order 
to answer the question, we use the Ornstein-Uhlenbeck 
process (OU process) as a simplest way to generate cor- 
related random variables with the correlation length L 
and check whether or not the coarse-grained quantities 
like Eqs.([T]) and @ follow lognormal distributions in the 
range R/L ~ 1. Our numerical results suggest that the 
answer to the question is yes. Then we study, by analyt- 
ical calculations of moments, further details on how the 
OU data become close to the lognormal variables. 

We believe that this simple approach using the stochas- 
tic process has some value since reproducing the large- 
scale lognormality by direct numerical simulations of the 
Navier-Stokes equations can be computationally quite ex- 
pensive. It requires a very large domain size as compared 
with the correlation length. 

In statistics, taking logarithm of a positive random 
variable is known as a common way of symmetrizing 
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transformation, which makes the skewness of the trans- 
formed variable closer to zero In this language, the 
large-scale lognormality suggests that, even with correla- 
tion, the log transformation can produce near Gaussian 
behaviour successfully already when the averaging scale 
R is of the order of the correlation scale L. The present 
study can be interpreted as a model study how the log 
transformation begins to work for the correlated random 
variable by changing the coarse-graining scale. 

This paper's organization is the following. In Sec. [TTJ 
we present numerical data of the OU process reproducing 
the large-scale lognormality. We then in Sec. IIIII provide 
analytical calculations of moments of the OU process to 
study how they are close to those of the lognormal dis- 
tribution. We provide summary and discussion in Sec. 

El 



II. NUMERICAL EXPERIMENT OF THE 
ORNSTEIN-UHLENBECK PROCESS 

We begin here by listing basic properties of the OU pro- 
cess (see, e.g., @ and also [§] in the context of turbulence 
modeling). The process is described by the Langevin 
equation 



dX(t) 1 . . , — , . 



(3) 



where the Langevin noise £(t) is a Gaussian white noise 
having the ensemble-averaged mean and variance 



<£(*)> = o, mm) = s(t-t'). 



(4) 



The linear addition of the Gaussian uncorrelated noise 
£(t) makes the solution X(t) to Eq.Q to be a Gaus- 
sian random variable. However the first term of the right 
hand side brings temporal correlation. Namely, by writ- 
ing the initial position as xq, X(t) is characterized as 
a correlated-in-time Gaussian random variable with the 
mean and variance 

(X(t)) = xoe-'/^, (5) 
({X(t)-(X(t))} 2 ) = «T L (l-e- 2 ^). (6) 

The correlation function can be calculated analytically 
as 

{lX(t)~(X(t))}{X(t + S )-(X(t + S ))}) = 

([X(t)-(X(t))] 2 )c-^ T \ (7) 

Therefore the integral scale of the OU process is given as 

[(X(t)-(X(t))}[X(t + S )-(X(t + S ))]) 



([x(t)-(x(tw 



-ds = T L .(8) 



If the large-scale lognormality observed in turbulence 
| is a property of correlated random variables, it is then 
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FIG. 1: (Color online) Probability density functions (PDF) 
of (a) X\ (Eq.©) and (b) its logarithm Xa.X\. The 
random variables are normalized to have zero mean and 
unit standard deviation. The averaging scales here are 
R = IOOOTl , IOOTl , IOTl , T L , Tl/ 10, which correspond to the 
curves from top to bottom. The PDFs are shifted vertically 
by being multiplied by factor 0.25 for clarity. The solid curve 
denotes the Gaussian distribution. 
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FIG. 2: (Color online) Same as Fig. [T] but for the 
squared increment (a) 8X^ R (Eq. (|10[0 and (b) its loga- 
rithm lnX^. Here r = T L /100. The averaging scales are 
R = IOOOTl , IOOTl , IOTl , T L , Tl/ 10, which correspond to the 
curves from top to bottom. 



likely that the following coarse-grained data of the OU 
process 



t+R/2 



t-R/2 



X(t) 2 dt, 



(9) 



I rt+R/2 

1 "> - ' [X(t + r)-X{t)] 2 dt,(lO) 



sx^ R (t) = 



R-r 



t-R/2 



are lognormally distributed when the averaging scale R 
is in the similar range, such as 1 < R/Tl < 100. Equa- 
tions © and (TTU1) correspond to the turbulence quantities 
Eqs.([T]) and ((2J, respectively. 

We now check whether or not the random variables 
(IH) and (fT0|) follow lognormal distributions by doing nu- 
merical simulation of the OU process (O- We use the 
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numerical method called exact updating formula of the 
OU process proposed in Ref.[10| with the parameters 
x = 0.0, Tl = 1.0, k = 0.50, At = 1.0 x 10" 3 . The 
integrals in Eqs. © and (ITU1) are numerically calculated 
as 



Nr/2 



Nr ^ 



X(t 



n+j/ 



(11) 



10 



-—N R /2 
Nr/2 



j=-N R /2 



15 
10 



(12) 



where N R = R/At,N r = r/At and t k = kAt. Fig- 
ures Q] and [2] show probability density functions (PDF) 
of X R and 8X 2 R with or without taking logarithm for 
various averaging scales R. The number of samples for 
the R = 1000T L case is 1.34 x 10 6 (for smaller R cases, 
the number is much larger). For the largest values of 
R = 1000Ti shown here, the PDFs are close to Gaussian 
distributions as a consequence of the central limit theo- 
rem. As we decrease the averaging scale R to the correla- 
tion scale Tl, the distributions of X R and SX 2 R deviate 
from the Gaussian distribution. In contrast, the log vari- 
ables lnX^ and \nSX 2 R remain to be nearly Gaussian 

as shown in FigsQJb) and[2jb). Hence X R and SX 2 R are 
lognormally distributed in this range of R. This behavior 
is similar to the turbulence data analyzed in Ref. [6|. In 
addition, for the squared increments SX 2 R with various 
r < Tl, qualitatively same results are obtained. 

The approach to the lognormal distribution can be ob- 
served more quantitatively by looking at how the skew- 
nesses and flatnesses of the log variables \nX R ,\nSX 2 R 
change as functions of R (the skewness S(z) and flatness 
F{z) of a random variable z are defined as S(z) = ((z — 
(z)) 3 )/[V(z)} 3 / 2 , F{z) = ({z- (z)) 4 )/[V(z)} 2 , where V(z) 
is the variance ((z — (z)) 2 )). In Fig. [31 it is seen that the 
moments of the variables with logarithm \nX R , \n5X 2 R 
approach already around R = Tl to the values of the 
Gaussian distributions (5, F) = (0, 3) whereas the mo- 
ments of the variables without logarithm are still different 
from them. This is consistent with the behavior observed 
in Figs H] and H 

For the increments SX 2 R , the skewness and the flat- 
ness can depend on the difference r. Indeed a clear in- 
dependence is seen in FigJU However, the fast conver- 
gence to the Gaussian of the \n5X 2 R around R ~ L is 
not affected by this dependence. In fact, these graphs of 
the different rs can be collapsed to one curve by normal- 
izing R with a different correlation scale from Tl as we 
will show at the end of the next section. 

In summary, we observe that the large-scale lognor- 
mality holds also for the OU process as in the turbulence 
case Q . Further details on the behavior of the skewness 
and flatness of the OU process are studied analytically 
in the next section. 
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FIG. 3: (Color online) The skewness of X R (a) and of SX^ R 
(r = Tl/100) (b); the flatness of X% (c) and of 5X* R (d) with 
or without logarithm. They are plotted versus the averaging 
scale R. The horizontal lines correspond to the skewness and 
the flatness of the Gaussian distributions. The thick curves 
are asymptotic expressions of the moments of X R and SX^ R , 
Eqs. (|25H25)) and (|28H29p . The dashed curves correspond the 
exact moment expressions (not involving any asymptotic ar- 
gument). 
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FIG. 4: (Color online) The skewness (a) and flatness (b) of 

ements SX^ R and \n5X; 
of r: r = Tt/IQ, T L /50, Tl/100. 



the increments 5X^ R and \nSX^ R for three different values 



III. ANALYTICAL EXPRESSIONS OF THE 
MOMENTS OF THE ORNSTEIN-UHLENBECK 
PROCESS 



The OU process is amenable to analytical calculation 
in many ways. However it is difficult to determine an- 
alytically the PDF's of the quantities of our interest, 
Eqs.© and (fTU)) , or all their moments, since they in- 
volve integrations. Instead we calculate the low-order 
moments, namely the skewness and the flatness of X(t) 2 R 
and SX(t) 2 R . For them, full analytic results are obtained 
but we present only the asymptotic expressions of them 
for R/L 3> 1 since they are sufficient for our purpose. 

Then we compare them with the expressions of the 
skewness and flatness of the lognormal distribution with 
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the mean e and the variance v. 



expressions: 



Sln = Vp(3 + P ), (13) 
-Fln = 3 + p(16 + 15p + 6p 2 + p 3 ), (14) 

where p — v/e 2 . As we have seen, the PDF's of the 
variables X(t) 2 R , SX(t) 2 R become close to the lognormal 
distributions for R/Tl ~ 1. We study this behavior by 
comparing the skewness and flatness of X(t) 2 R , 5X(t) 2 R 
with Sln and Fln- However, as we shall see below, it 
turns out that small corrections exist, which may be hard 
to be detected in numerical calculations in the previous 
section. 

In our analytical calculation of the moments, we first 
put the formal solution of the Langevin equation (|3|) 



X(t) =e- t / TL x a + 



[ e-^- s ^ TL VM(s)ds. (15) 
Jo 



into the coarse-grained variables Eqs.© and (flO]) . For 
the correlation or the variance, we take a suitable powers 
of them and take the ensemble average with respect to the 
Langevin noise by using Eq.Q. So that we can reduce 
their calculations into multiple integrals of certain expo- 
nential functions. For the skeweness and the flatness, we 
use more compact integral representations proposed in 
Ref.pj]], which is explained in Appendix. To calculate 
the resultant integrals we use a symbolic computation 
software, Maple. The integrated result still contains de- 
pendence on time t. For example, the mean of X(t) R 
reads 



(X(t) 2 R ) = nT L 



2ir 



)e 



-2t/T L 



(16) 



Remember that this time t is much larger than the cor- 
relation time Tl to be consistent to the simulations in 
Sec|n] We hence regard <&~ 1 I Tl = 0. From now on, we 
drop the dependency on t of all the moments of X(t) R 
and 8X(t) 2 R calculated here. This does not affect our 
study since the OU process becomes steady for t >• Tl- 
We now list the expressions of the moments. Here we 
write A = R/Tl for brevity. For X R , the mean, variance, 
skewness, and flatness are respectively, 



E{X R ) 
V(X%) 

S(X%) 
F(X 2 R ) 



(kT l ) 2 [2A- 1 



+ A- 2 (e 



-2A 



(17) 
(18) 



= 12V' [e -l+A(e" JA + l)], (19) 



= 3 



6WV 2 [c- 4A 



28e _2A - 29 



-16AV 2A 



40Ae" 2A + 20A], (20) 



where Wq = 2A 1 +A 2 (c 2A — 1). For the coarse grained 
increment SX 2 R , by writing r/Tj, = A, we give the results 
in the leading order of A -1 = (R/Tl) -1 to avoid lengthy 



E{5X 2 R ) = 2 K T L e- x (e x - 1), 



(21) 



V{8Xl R ) 



(2KT L ) 2 A- 1 e- 2A [3e 2A - 4(A + l)c A 

+2A + 1] , (22) 



S(SX, 



r,R) 



3A- 2 e" 3A [l«e 3A - 5(A 2 + 3A + 3)e 2A 

+2(4A 2 + 6A + 3)e A - (3A 2 + 3A + 1)] , (23) 
F(6X 2 R ) = 3 + A^ 3 e^ 4A W/f 2 [525e 4A 

-(56A 3 + 336A 2 + 840A + 840)e 3A 
+(224A 3 + 672A 2 + 840A + 420)e 2A 
-(216A 3 + 432A 2 + 360A + 120)e A 
+ (64A 3 + 96A 2 + 60A+15)], (24) 

where Wi = A _1 e _2A [3e 2A - 4(A + l)e A + 2A+ 1]. Notice 
that we obtain analytic expressions of S and F, which 
are shown as the dashed curves in Fig|3] 

With these asymptotic expressions, we next re-write 
the skewness S and flatness F as functions of the variance 
over the squared mean p — V/E 2 to compare them with 
the lognormal ones Sln and Fln- 

For X R , in the limits of A = R/Tl — > oo, we have 



S(X 2 R ) ~ + 4~p 3 



where 



16 P 



Ft ) ' 1 - /> ( 15 + |p - |p 2 



V(X 2 R ) 2T L (T L 



E{X 2 ) 2 



R 



R 



(25) 
(26) 

(27) 



The asymptotic expression Eqs. (l2~5|) and (|2"6"|) agree with 
data shown in FigE]for R > VlQT L . 

For the velocity increments SX 2 Rl we have p = 
V/E 2 = A- 1 [3e 2A - 4(A + l)c A + 2A + l]/(e A - l) 2 . By 
taking the limit of A = r/T^ — > 0, we obtain expressions 



!)() 



S(SX^ R ) ~ o + vpx - 



F(SX 2 R ) ~ 3 + px 



40' 
1359 

ISo"' 



where 



P : 



E(5X 2 M Y 



4A _ 4r 
3A ~ 3R' 



(28) 
(29) 

(30) 



The equations (l28t and (|29|) agree with the data shown in 
Fig|3]for R > 10 _3 / 2 Tl (covering all data points), thanks 
to the small value r = 10 _2 Fl. 

We now compare the analytic expressions of the OU 
quantities Eas. (|25H26)l and (p8] - [29f with the log-normal 
ones Eqs. p^HT^|) . We can observe the following: (i) In- 
deed as R — !• oo (p -t 0), S and F of both X R and SX 2 R 
tend to the values of Gaussian distribution, S = and 
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F = 3. Before reaching this state, we see an approach 
to the lognormal distribution, (ii) The sub-leading terms 
of S and F have the same powers of p as Sln and -Fln ; 
which is in favor of the large-scale lognormality of the OU^ 
process. However, the constant in the subleading term 
is slightly different from those of the lognormal distri- 
bution. In this sense, the coarse-grained quantities X R 
and 5X 2 R of the OU process becomes nearly lognormally 
distributed when R/Tl is large but not exactly so. 

So far we focus on how the moments vary as a func- 
tion of R/Tl, where R is normalized by the correlation 
scale Tl of the OU process. We close this section with 
a digression by pointing out another normalization scale 
of the problem. This is indicated by the above analytic 
results. For X R) let us go back to the definition Eq.©. 
The integral scale of the integrand of Eq. ([9]) is 



-(x 2 ) = 



; oo ([x 2 (<)-(x 2 )][x 2 (t + s )-(x 2 )])d S 



([x 2 (t)-(x 2 )] 2 ) 



Tl 
2 ' 



(31) 



For SX 2 R , the integral scale of its integrand defined sim- 
ilarly is calculated as 



t(6X 2 ) = T L 



3c 



2A 



4(A + l)e A + 2A + l 
4(e A - l) 2 



(32) 



for small A = t/Tl- Notice that the p variables defined 
in Eqs. ((27|) and (|30| contains the corresponding integral 
scales of the integrands (|3"T|) and ([32|) . namely, p = At/R 
for both X R and 6X* R . Hence it implies that a suit- 
able normalization scale of R for collapsing the moment 
data plotted in Figs H] and |U are respectively the integral 
scales t(X 2 ) and t(SX 2 ) of the integrands. Indeed this 
normalization yields a better collapse as shown in FigjS] 
for both X R and 5X 2 R with different r. This suggests 
that, when we look into a coarse-grained quantity over 
scale R of a function fn{y) = (1/-R) Jy +R .fW)dy' of a 
correlated fluctuation y, the correlation scale (integral 
scale) of f(y) in certain cases can be a better character- 
istic scale than the correlation scale of the fluctuation y 
itself. We believe that this is the case not only for the 
OU process but also for a general fluctuation with a cor- 
relation. However this normalization blurs the fact that 
the large-scale lognormality of the OU process occurs at 
R/Tl ~ 1. This kind of normalization will be reported 
elsewhere. 
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FIG. 5: (Color online) The collapsed skewness (a) and flatness 
(b) of X\ and SX% R as functions of R normalized with t(X 2 ) 
and t(5X?) defined in ^STJ and (J32J) . 



and velocity increments become lognormally distributed 
fluctuations. We here anticipated that this large-scale 
lognormality is a property of correlated random variables 
in general. 

To test this idea, we took the Ornstein-Uhlenbeck pro- 
cess as a simplest means to generate correlated random 
variables. Our numerical simulation indicates that the 
OU process also exhibits the large-scale lognormality for 
about the same range of R/L as observed in turbulence. 
Then, to further investigate the approach to the lognor- 
mal distribution, we calculated analytically the skewness 
and flatness of the coarse-grained quantities of the OU 
process. The first subleading term of the asymptotic ex- 
pressions for large R/L of the moments revealed that 
indeed the moments behave nearly as those of the log- 
normal distribution. However there is small deviation 
from the lognormal distribution. 

Now we speculate on behavior of the higher order mo- 
ments than forth order. The analytical calculation of 
them for the coarse-grained OU variables becomes so 
complicated that we did not try even with Maple. For 
the lognormal variable ^ln j the moments can be written 
with p, the ratio of the variance to the squared mean, as 



(q) _ ([gLN - (^Ln)] 9 ) 



LN 



[U(x LN )K 2 

9 



k=0 



(P+1) 



jfc(fc-l) 



(-1) 



q—k 



(33) 



where (?) is the binomial coefficients. By expanding 
Eq. (|33[> with p, the moment can be written as 



(?) 

LN 







+ p i«(9-2) ( g; dd), 



(q - I)!! + f[ q) p + f^'p 2 + ... + pUi-V (q: even) 



(34) 



IV. SUMMARY AND DISCUSSION 

This study was motivated by the large-scale lognor- 
mality of turbulence that was recently observed experi- 
mentally in grid, boundary- layer, and jet turbulences Q. 
In this lognormality, the correlation scale plays a pivotal 
role. Namely when the averaging scale to the correlation 
scale are of order unity, the averaged squared velocity 



Here the terms of the negative powers of p are shown to 
vanish and s±,f[,... are suitable coefficients, which 
can be written with the binomial coefficients. The p° 
terms in Eq. flMl) . namely and (q — I)!! = (q — l)(q — 
3) • • • 5 • 3 • 1, correspond to the values of the Gaussian 
distribution. We speculate that, for the coarse-grained 
quantities of the OU variable, their higher-order mo- 
ments have the similar expression as Eq. (|34[) . As we saw 
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in the previous section, the behavior for small p (corre- 
sponding to R/Tl being larger than unity) is relevant for 
the large-scale lognormality. Their first term should co- 
incide with that of Eq. ([M|) because of the usual central 
limit theorem. It is tempting to speculate also that the 
second term has the same power of p but the value of the 
the coefficient, Sj 9 or f[ q ^ are different. We have no idea 
how much they are different. 

In relation to the real turbulence data, the OU pro- 
cess is not a good representation as a whole since, for 
example, it does not have deviation from being Gaussian 
(intermittency effect) and equivalence of the energy cas- 
cade nor the energy dissipation rate. However, as long 
as we focus on the large-scale fluctuations of turbulence 
where the single point velocity or the velocity increments 
are close to being Gaussian, the OU process serves a use- 
ful and analytically tractable model. In this study we 
regarded that the correlation at the integral scale is the 
most important aspect to be modeled by the OU process. 

Here we have seen that the correlation plays the es- 
sential role in the near lognormal behavior of the coarse- 
grained positive quantities Eqs.© and (|T0|. Roughly 
speaking, this near lognormality around the correlation 
scale may be regarded as an intermediate state, or "a 
rule of thumb" before the central limit theorem holds 
with much larger averaging scale R. For the turbu- 
lence cases studied in [6|], similar mechanism is likely to 
work. Other examples of lognormal behavior involving a 
coarse-graining av erag e include cosmological density fluc- 
tuations (see e.g., pjj])) which may have a similar struc- 
ture as studied here. 



V. APPENDIX 

Here we explain briefly how we calculate the moments 
of the coarse-grained quantities. In particular we use the 
Rice's method given in Section 3.9 of Ref.fliT] to reduce 
the number of the multiple integrals. For illustration of 
the method, let us take as an example the second-order 
moment of X R (t) 



t + 4 



dt 1 



t-R/2 



t+4 



t-R/2 



dt 2 (X^X'fo)). 



(35) 



The idea of the Rice's method is to write the mo- 
ment (X 2 (ti)X 2 (t2)} in terms of the correlation function 
(X(*i)X(i 2 )) = tp(t 2 - ii) by noticing that X(t t ) and 
X(t 2 ) are multivariate Gaussian variables whose covari- 
ace matrix is known completely In principle, any mo- 
ment of X(t±), X(t 2 ), ■ ■ ■ , X(t n ) can be expressed with 
ip(ti — tj). Specifically, with the correlation function for 
large t 



^{t) = (X(t)X(t + T)) = K T L e 



\t\/Tl 



(36) 



the variance and other central moments can be written 
as follows: 



t+4 



(ft, 



t+4 



dt 2 ip 2 (ti-t 2 ) 



R 2 Jo 
K 2 Tl r 

4i? 2 



{R- x)ip 2 (x)dx 



2T L R + Tl{e 



-2R/T L 



1) 

(37) 



O pt-\--% ft-\--2 rt-\--2 

([X 2 R (t) - (X 2 )} 3 } = — I dtl I dt 2 I dt 3 

Jt — -j- Jt — -j- Jt — 

xV(*i-*aW*a-*3M*3-*i) 
48 f R 

= -53 / dx{R-x)ip{R-x) 
R Jo 

x f dy^(x-y)^(y) (38) 
Jo 

([X 2 R {t)-(X 2 R )f) = 3((X 2 R (t)-(X 2 R )) 2 ) 2 

48 r t+ -2 r t+ - r t+ - r t+ - 

+SI / dt i / dt 2 / dt 3 / dU 
R ./* « Ji-S, Jt-& Jt.-S. 



xip(t 2 - ii)V>(*3 - *i)V(*4 - h)4>(U - t 3 ) 
3((X 2 (t) (X 2 R )) 2 ) 2 



pR fR 

/ dx(R - x) I dynp(y)ip{x - y) 
Jo Jo 
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x / dzi\)(x - z)tJj(z). (39) 
Jo 

Here we change variables in the integrals by using the 
symmetry of ip to reduce the double integral to a single 
integral (for details, see [ill]). The integrals (|3"5)) and (|59")l 
are calculated analytically with the software Maple. 

Concerning the increments, we write its correlation 
function for large t 



i/j r ( T ) = (S r X{t + T)S r X{t)) 



= kT l 



,-\r\/T L{ 2_ e -r/T L) 



\/Tl 



(40) 



Using the same argument as in the case of X R (t), we can 
express the central moments with the correlation ip r as 
follows: 



([SX 2 R (t) - (SX 2 R (t)}} 2 ) 



([SX 2 (t) (SX'it))} 3 ) = 



R-r 



(R r) 2 J 

xtp 2 (x)dx 



(R-r- x) 



/■-R-r 

/ dx(R — r — x) 
Jo 



(R - rf 

X / dytp r (x - y)lpr(x)lpr{y) 

Jo 

(41) 
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([6X* R (t) - (SX^ R (t))} 4 ) Acknowledgments 

= 3{{SX* R (t) {5Xl R {t))] 2 ) 2 
768 rR ~ 



(R 

R—r py 



r-tt—r 

/ dx(R — r — x) w e gratefully acknowledge discussion with H. Mouri. 
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